IIIIIIM 

US007701561B2 

(i2) United States Patent (io) Patent no.: us 7,701,561 B2 

Zou et al. (45) Date of Patent: Apr. 20, 2010 


(54) POLARIZATION IMAGING APPARATUS 

(75) Inventors: Yingyin Kevin Zou, Lexington, MA 
(US); Qiushui Chen, Andover, MA 
(US); Hongzhi Zhao, Malden, MA (US) 

(73) Assignee: Boston Applied Technologies 

Incorporated, Woburn, MA (US) 

( * ) Notice: Subject to any disclaimer, the term of this 

patent is extended or adjusted under 35 
U.S.C. 154(b) by 507 days. 

(21) Appl.No.: 11/701,738 

(22) Filed: Feb. 2, 2007 

(65) Prior Publication Data 

US 2007/0296958 Al Dec. 27, 2007 

Related U.S. Application Data 

(60) Provisional application No. 60/772,779, filed on Feb. 

13, 2006. 

(51) Int.Cl. 

G01N 21/00 (2006.01) 


(52) U.S. Cl 356/73 

(58) Field of Classification Search 356/72, 


356/73, 322 

See application file for complete search history. 

(56) References Cited 

U.S. PATENT DOCUMENTS 

5,412,500 A * 5/1995 Fergason 349/76 

5,847,394 A 12/1998 Alfanoetal. 

6,437,856 B1 8/2002 Jacques 

6,700,694 B2 3/2004 Zouetal. 

6,746,618 B2 6/2004 Lietal. 


6,798,514 B2 9/2004 Daniels 

6,816,261 B2 11/2004 Patel etal. 

6,890,874 B1 5/2005 Lietal. 

2003/0067667 Al * 4/2003 Wang etal 359/323 

OTHER PUBLICATIONS 

S. Rahmann et al., “Reconstruction of specular surfaces using polar- 
ization imaging,” IEEE Computer Society Conference on Computer 
Vision and Pattern Recognition, vol. 1, pp. 149-155, 2001. 
S.G.Demos, et al., “Deep subsurface imaging in tissues using spectral 
and polarization filtering,” Optics Express, vol. 7, No. 1, pp. 23-28, 
2000. 

H. Zhao, et al., “Polarization imaging sensor for cell and tissue 
imaging and diagnostics,” Proc. of SPIE vol. 6380, pp. 638006-1-11, 
2006. 

* cited by examiner 

Primary Examiner — Tarifur Chowdhury 
Assistant Examiner — Abdul lahi Nur 

(57) ABSTRACT 

A polarization imaging apparatus measures the Stokes image 
of a sample. The apparatus consists of an optical lens set 11, 
a linear polarizer 14 with its optical axis 18, a first variable 
phase retarder 12 with its optical axis 16 aligned 22.5° to axis 
18, a second variable phase retarder 13 with its optical axis 17 
aligned 45° to axis 18, a imaging sensor 15 for sensing the 
intensity images of the sample, a controller 101 and a com- 
puter 102. Two variable phase retarders 12 and 13 were con- 
trolled independently by a computer 102 through a controller 
unit 101 which generates a sequential of voltages to control 
the phase retardations of VPRs 12 and 13. A set of four 
intensity images, I 0 , I l5 1 2 and I 3 of the sample were captured 
by imaging sensor 15 when the phase retardations of VPRs 12 
and 13 were set at (0,0), (jt,0), (jt,jt) and (jc/2, jc), respectively 
Then four Stokes components of a Stokes image, S 0 , S l3 S 2 
and S 3 were calculated using the four intensity images. 

10 Claims, 5 Drawing Sheets 
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POLARIZATION IMAGING APPARATUS 

CROSS REFERENCE TO RELATED 
APPLICATION 

5 

This application claims the benefit of Provisional Patent 
Application Ser. No. 60/772,779 filed Feb. 13, 2006. 

STATEMENT OF GOVERNMENTAL INTEREST 

to 

This invention was made partially with Government sup- 
ports under grant no NNJ05JC13C awarded by NASA. The 
government has certain rights in this invention. 

BACKGROUND OF THE INVENTION 15 

1 . Technical Field 

The invention relates to optical imaging apparatus, particu- 
larly related to polarization imaging apparatus, more specifi- 
cally related to a high speed polarization imaging system that 2 o 
can measure Stokes components from a sample at different 
optical wavelengths. 

2. Technical Background 

The light scattered by a tissue has interacted with the ultra- 
structure of the tissue, which imprinted some intrinsic prop- 25 
erties of the tissue. Tissue ultra structure extends from mem- 
branes to membrane aggregates to collagen fibers to nuclei to 
cells. Photons are most strongly scattered by those structures 
whose size matches the photon wavelength. It has been dem- 
onstrated that light scattering can provide structural and func- 30 
tional information about the tissue. One important biomedi- 
cal application of optical imaging and spectroscopy is non- 
invasive or minimally invasive detection of pre-cancerous and 
early cancerous changes in human epithelium, such as dys- 
plasia or carcinoma in situ. 35 

In recent years there has been an increasing interest in the 
propagation of polarized light in randomly scattering media. 
The investigation of backscattered light is of particular inter- 
est since most medical applications aimed at the in-vivo char- 
acterization of biological tissue rely on backscattered light. 40 
By recording the spatially dependent response of a medium to 
a polarized point source, one may obtain information about 
the scattering particles that are not accessible to mere inten- 
sity measurements. A diagnostic imaging modality based on 
near-infrared (NIR) radiation offers several potential advan- 45 
tages over existing radiological techniques. First, the radia- 
tion is non-ionizing, and therefore reasonable doses can be 
repeatedly employed without harm to the patient. Second, 
optical methods offer the potential to differentiate between 
soft tissues, due to their different absorption or scatter at NIR 50 
wavelengths that are indistinguishable using other modali- 
ties. And third, specific absorption by natural chromophores 
(such as oxy -haemoglobin) allows functional information to 
be obtained. NIR imaging research has focused on a variety of 
possible clinical applications. 55 

The fact that the polarization state of the light contains 
useful information has been shown in many literatures, for 
example. Rahmann and Canterakis in “Reconstruction of 
specular surfaces using polarization imaging,” IEEE Com- 
puter Society Conference on Computer Vision and Pattern 60 
Recognition , vol. 1, pp. 149-155, 2001, describe how the 
polarization state of light can be used for the reconstruction of 
specular surfaces to determine the shape of 3-D objects. They 
use the fact that light reflected by dielectrics and metals 
becomes partially linearly polarized and that the direction of 65 
polarization depends on the orientation of the reflecting sur- 
face. Demos and Alfano in “Deep subsurface imaging in 
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tissues using spectral and polarization filtering,” Optics 
Express , vol. 7, no. 1, pp. 23-28, 2000, demonstrate a tech- 
nique based on polarization imaging that allows for optical 
imaging of the surface as well as structures beneath the sur- 
face. 

There have been several polarization imaging schemes 
published. InU.S. Pat. No. 6,437,856, inventor Steven Louis 
Jacques, which is sued Aug. 20, 2002, andU.S. Pat. No. 5,847, 
394, inventors Robert R. Alfano, et al., which issued Dec. 8 , 
1998, a set of measurements at different polarization orien- 
tations are taken to render a new image that is independent of 
the light reflected from the surface of a tissue sample and that 
is dependent of the light scattered from deep tissue layers. 
Especially, a linearly parallel polarized light is used for illu- 
mination and two images are acquired, one image selecting 
linearly parallel (Par) polarized light (i.e., parallel to the light 
source-tissue-camera plane) and one image selecting linearly 
perpendicular (Per) polarized light (i.e., perpendicular to the 
light source-tissue-camera plane). A new image is obtained 
by (Par-Per) or (Par-Per)/(Par+Per). 

Many polarimetric sensing technologies have been devel- 
oped to capture the Stokes polarization information. People 
use rotating retarder, rotating polarizer to obtain Stokes 
parameters from several successive frames of image. In these 
mechanically-rotation approaches, however, the exact spatial 
registration between various frames is difficult due to changes 
of viewing angle and the polarization properties during the 
period of between successive frames in most target detection 
applications. 

In U.S. Pat. No. 6,798,514, inventor James Maurice 
Daniels, which issued Sep. 28, 2004, the light passes through 
two liquid crystal waveplates and a polarizing filter before 
falling on a light sensitive element device to measure intensity 
of light in four different polarizations. However, the same 
difficulty also applied to the approach using liquid crystal for 
electro-optic polarization modulation, which usually has 
response time in the order of 100 ms. Another drawback for 
the liquid crystal approach is that polarization-rotating liquid 
crystals are not commercially available in many important IR 
bands. 

SUMMARY OF THE INVENTION 

It is therefore the main objects of the invention to provide 
a polarization imaging apparatus which can consistently pro- 
vide high speed and broadband Stokes images in real-time. 

The first key feature of the polarization imaging apparatus 
is that it provides fast adjustable optical phase control through 
high speed electro -optical variable phase retarders made from 
transparent electro-optic ceramics. 

The second key feature of the invention is that the polar- 
ization imaging apparatus can be worked at a wavelength 
from 400 nm to 7 pm, using electro-optic ceramics. 

The third key feature of the invention is that a simplified 
time sequential control and measurement algorithm which 
provides a fast measurement and converts sequential intensity 
images to Stokes images. 

A preferred embodiment of the polarization imaging appa- 
ratus is shown in FIG. 1, which consists of a polarization 
imaging sensing unit 100, an electronic driver 101 and a 
computer 102 . The imaging sensing unit 100 consists of an 
optical lens set 11 , a linear polarizer 14 with its optical axis 
18 , a first variable phase retarder (VPR) 12 with its optical 
axis 16 aligned 22.5° to axis 18 , a second variable phase 
retarder 13 with its optical axis 17 aligned 45° to axis 18 , and 
a imaging sensor 15 for sensing the intensity images of the 
sample. The imaging sensor 15 can be a focal plan array 
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(FPA) detector, a CMOS or a CCD sensor. Two VPRs 12 and 
13 were controlled independently by a computer 102 through 
a controller unit 101 which generates a sequential of voltages 
to control the phase retardations of VPRs 12 and 13. A set of 
four intensity images, I 0 , I 1? I 2 and I 3 of the sample were 5 
captured by imaging sensor 15 when the phase retardations of 
VPRs 12 and 13 were set at (0,0), (jt,0), (jt,jt) and (ji/2,jt), 
respectively Then four Stokes components of a Stokes image, 

S 0 , S l5 S 2 and S 3 can be calculated by computer 102 using the 
following formula: 1C 


' So = 0.5x(/ 0 + /i) (1) 

S { =0.5x(/ 0 -/i) 

' S 2 = 0.5x(/ 2 -S 0 ) 15 

,S 3 =0.5x(/ 3 -So) 

In one aspect, the present invention includes the ferro- 
electric materials, such as PLZT, La modified PMN-PT, and 20 
La modified PZN-PT to form the electric voltage controlled 
ferro-electric variable phase retarder. 

Referring to FIG. 5 is another preferred embodiment of the 
polarization imaging apparatus consists of an illumination 
unit 50, a polarization imaging sensing unit 500, a controller 
unit 101 and a computer 102. The said illumination unit and 
imaging sensing unit are aligned in a reflection configuration. 

In addition to a light source, the illuminator may consist of 
parts such as a collimator, a polarizer or polarization modu- 
lator, and a fixed wavelength filter or tunable wavelength 
filter. Tunable optical band-pass filter is used to produce 
desired wavelength for the imaging, from visible to infrared 
(IR). The illuminator can produce a linear or arbitrary polar- 
ized light beam. 

The illumination unit 50, and a polarization imaging sens- 
ing unit 500 in FIG. 5 can be further aligned in a transmission 
configuration, as illustrated in FIG. 6. 

Additional features and advantages of the invention will be 
set forth in the detailed description which follows, and in part 40 
will be readily apparent to those skilled in the art from that 
description or recognized by practicing the invention as 
described herein, including the detailed description which 
follows, the claims, as well as the appended drawings. 

It is to be understood that both the foregoing general 45 
description and the following detailed description present 
embodiments of the invention, and are intended to provide an 
overview or framework for understanding the nature and 
character of the invention as it is claimed. The accompanying 
drawings are included to provide a further understanding of 50 
the invention, and are incorporated into and constitute a part 
of this specification. The drawings illustrate various embodi- 
ments of the invention and together with the description serve 
to explain the principles and operations of the invention. 

55 

BRIEF DESCRIPTION OF THE DRAWINGS 
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FIG. 5 is a schematic view of a preferred embodiment of 
reflective polarization imaging apparatus in accordance with 
the present invention; 

FIG. 6 is a schematic view of a preferred embodiment of 
transmissive polarization imaging apparatus in accordance 
with the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED 
EMBODIMENT 

It is to be understood that the invention is not limited in its 
application to the details of the particular arrangement shown 
since the invention is capable of other embodiments. Also, the 
terminology used herein is for the purpose of description and 
not of limitation. 

FIG. 1 shows a preferred embodiment of the polarization 
imaging apparatus in accordance with the present invention. 
It consists of an imaging sensing unit 100, an electronic driver 

101 and a computer 102. The imaging sensing unit 100 fur- 
ther consists of an optical lens set 11, a polarizer 14 with its 
optical axis 18, a first variable phase retarder (VPR) 12 with 
its optical axis 16 aligned 22 . 5 ° to axis 18, a second variable 
phase retarder 13 with its optical axis 17 aligned 45 ° to axis 
18, a imaging sensor 15 for sensing the intensity images of the 
sample. The imaging sensor 15 can be a focal plan array 
(FPA) detector, a CMOS or a CCD sensor. The imaging 
sensing unit then is assembled inside a metal host (not 
shown). Two VPRs 12 and 13 were controlled independently 
by a computer 102 through a controller 101 which generates 
a sequential of voltages to control the phase retardations of 
the VPRs 12 and 13. The intensity images were captured 
through the imaging sensor 15, and processed by computer 

102 to convert the sequential intensity images to Stokes 
images. 

Using Mueller’s calculus, we can obtain the Stokes param- 
eters of an input light 10 by measuring the output light inten- 
sity of the sensor 15. The transformation done by the phase 
retarders can be described in the following equation, 

Sour(So,ou^l,ou^2,ou^3,ou t ) T =K 2 *M r i (2) 

where S IW =(S 0 S : S 2 S 3 ) r is the unknown input Stokes vec- 
tor of the target 10, which consists of four Stokes compo- 
nents. The parameter S 0 is the intensity of the light (including 
un-polarized component), S : is the intensity difference 
between horizontal linearly polarized and the vertical linearly 
polarized components, S 2 is the intensity difference between 
45 ° linearly polarized and the 135 ° linearly polarized com- 
ponents, and S 3 is the difference between left- and right- 
circularly polarized intensities. 

The quantities M rl , and M r2 are the Mueller matrices rep- 
resenting two variable phase retarders 12 and 13, respectively. 
These are given for the specific orientation by 


FIG. 1 is the schematic of a preferred embodiment of 
polarization imaging apparatus with variable retarders in 
accordance with the present invention; 

FIG. 2 shows the measurement steps for polarization imag- 
ing apparatus in accordance with the present invention; 

FIG. 3 shows the schematic drawing of a high speed vari- 
able retarder using a electro-optic ceramic; 

FIG. 4 shows measured responses of a high speed variable 
retarder; 
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where 61 and 62 are the two phase retardance introduced 
by the two variable phase retarders 12 and 13, respectively. 10 
To be able to measure the four components of the input 
Stokes vector, we need to set the variable phase retarders to at 
least four different states labeled as k. In the following, we 
detail the linkage between input state Stokes parameters and 
measured light intensity. The Mueller matrix for the cascaded 1 5 
variable retarders can be generally cast into an alternative 
form so that the state of polarization incident on the polarizer 
can be expressed as, 
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The actually detected signal through the polarizer is essen- 
tially 


1 

^ — 2 i^0,out + 


30 


Without worrying about a global factor, we throw away the 35 
factor of V 2 , thus the intensity can be expressed in terms of 
three of the Muller matrix elements, 


I-So,out+S 1 , out - m OoSo,in +m 1 A l,in +m 12^2, in +m 1 3$3,in ( 6 ) 

It is clear that the output intensity can be linked to the input 40 
Stokes parameters through generating four known states of 
polarization at the input so that we have, 
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(7) 45 


50 


where m tJ k is the Mueller matrix element at the k-th state of 
variable retarders, and I z is the intensity measured with the i-th 
state of variable retarder. From Equation (7), it can be noticed 
that once the matrix A is known, the input Stokes parameters 55 
can be found out the inverse the matrix so that, 


r So, in ' 
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The very question now is to find out the matrix A. It can be 65 
calculated from each Muller matrix in the system or, it can be 
done by generating four known states of polarization so that, 


where is the i-th Stokes parameter generate by the j -th 

retardance. Thus we are able to link the input state polariza- 
tion to the intensity measured at four retardance states. The 
Mueller matrix elements needed for finding the input state of 
polarization can be figured out through the use of four known 
state of polarization. 

The measurement procedure is based on specifying above 
general principle for practical applications. As the configura- 
tion shown in FIG. 1, we can find out the useful Muller matrix 
elements (according to equation ( 6 )) as follows: 

moo=l 

^cos 2 (81/2) -cos (82)-sin(8l)- sin (82)/V2 
m 12 =sin 2 (8l/2) 


w 1 3=-cos 2 (8l/2)-siii(82)-siii(8l)-cos , (82)/v^2 (10) 

where 61 and 62 are the retardance of two VPRs 12 and 13, 
respectively. Although their values could be many different 
combinations for finding out the Stokes vector, the following 
4-state procedure is considered to be the simplest yet most 
effective approach. 

Referring to FIG. 2 shows the procedure diagram to 
accomplish one measurement of Stokes parameters. 

Four intensity images (I 0 -I 3 ) are needed and taken sequen- 
tially, with two variable phase retarders 12 and 13 controlled 
by a sequence of voltages which were generated by controller 
1 01 . Phase retardation changes 21 and 22 of VPR 12 and VPR 
13, respectively. I 0 is taken with 0 volt applied to both 12 and 
13. 1 1 is taken with the half-wave voltage applied to 12 and 

0 volt applied to 13. I 2 is taken with is applied to both 
retarder 12 and 13. I 3 is taken with a quarter- wave voltage 
Vyt /2 applied to 12 and applied to 13. The four Stokes 
vector images S 0 -S 3 can be calculated from the intensity as 
follows: 


50 = 0.5x(/ 0 + /i) (11) 

51 =0.5x(/ 0 -/i) 

5*2 = 0.5 x (/ 2 - So) 

S 3 = 0.5 x (/ 3 -S 0 ) 


Table 1 listed the state of retardance of first and second 
variable phase retarders and the corresponding four intensi- 
ties measured in each equal length time slot. The intensities 
also can be expressed as a function of Stokes parameters 
S 0 -S 3 of the input Stokes vector according to equation (2). 
Table 1 also listed the formula for extracting the Stokes 
parameters S 0 -S 3 from these intensities. 

TABLE 1 


Time slot 1 2 3 4 

State of VPR1 (81) 0 % jr jt/2 

State ofVPR2 (82) 0 0 Jt jt 

Measured intensity I 0 h I 2 I 3 
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TABLE 1 -continued 


Time slot 

1 

2 

3 4 

Relation to Stokes 

Sq + Sj s. 

0,zn _S l,zn S w 

n + ^2 ,in ^0,in + 

Derived parameters 


' So = 0.5 X (I 0 + P) 

51 = 0.5 X (I 0 — Ii) 

5 2 = 0.5 x (I 2 - Sq) 
, S3 = 0.5 x (I 3 - Sq) 


With this simple converting formula, Stokes images can be 
obtained using very simple calculation from the intensity 
images captured by an imaging sensor. 15 

It should be pointed out that the switch of the relative 
positions ofVPR 12 and VPR 13 will not affect the outcome 
of the Stokes images according to equation (2). 

It is to be appreciated that the orientation configuration 
between the two variable phase plates 12 and 13 might be 20 
achieved with an intervening fixed waveplate, so that physical 
orientation of the two plates could be less constrained. Hence, 
the 22.5 degree orientation of the second variable phase 
retarder 13 with respect to the first variable phase retarder 12 
of claim 1 is achieved by any combination of relative physical 25 
rotation of the second variable phase retarder 13 with respect 
to the first retarder 12 and rotation of the polarization refer- 
ence frame between the two variable phase retarders 12 and 
13 by a fixed phase retarding plate (not shown but included as 
part of the present invention). 30 

Any other orientation configuration of first variable phase 
retarder 12 and second variable phase retarder 13 will also 
work, as long as the two variable retarders oriented at differ- 
ent angles, however may involve more difficult calculations. 

In addition to the simplified measurement algorithm, a key 35 
component in this polarization imaging apparatus is the vari- 
able phase retarder. The other means to change the phase 
retardation, rotating waveplates, Faraday rotators and Pock- 
el’s cell have been used in laboratories worldwide. These 
devices, however, are bulky, high power consumption, slow 40 
or expensive. 

A preferred variable phase retarder is made from an elec- 
tro-optic (EO) material. The general requirement for an EO 
variable phase retarder is that, when a voltage is applied, a 
predetermined phase shift is produced for the light beam. 45 
Also the material is isotropic with no voltage applied, there is 
no residual birefringence to cause the initial phase shift. Pref- 
erably, the material has a high electro -optic coefficient in 
order to reduce operating voltages to less than 500 volts. Also 
the mechanical characteristics allow formation of a bar or 50 
plate for use as the electrode. Of course, the material must be 
transparent at the wavelength of interest, e.g., between 450 
nm and 7 pm for most of the polarization imaging applica- 
tions. 

These requirements are satisfied by a class of ferroelectric 55 
complex oxides which 1) are optically isotropic; 2) have a 
Curie temperature less than about 490° C., so that electro- 
optic coefficients are high near room temperature; 3) have a 
diffusive phase transition, so that the temperature dependence 
of the electro-optic coefficients is lessened; and 4) which are 60 
not permanently poled by moderate electric fields since mate- 
rials with a low Curie temperature that become permanently 
poled are less stable. Example material systems include elec- 
tro-optic lanthanum modified PZT (PLZT), lead manganese 
niobate (PMN), a solid solution of lead zinc niobate and lead 65 
tantalate (PZN-PT), and a solid solution of lead manganese 
niobate and lead tantalate (PMN-PT). Besides PLZT, without 
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being an exhaustive list the following materials and their solid 
solution with or without lanthanum modification, may be 
used: Pb(Zr,Ti)0 3 , Pb(Mg,Nb)0 3 , PbTi0 3 , Pb(Zn,Nb)0 3 and 
PbTa0 3 . These materials are available from Coming Inc. 
(Corning, N.Y.) and Boston Applied Technologies, Inc. 
(Woburn, Mass.), referenced as OptoCeramic® materials. 
These materials are electrical controllable, solid-state, high 
speed, low insertion loss, and low fabrication cost. 

A preferred OptoCeramic® material, PMN-PT, for use in 
the present invention may be described by the general for- 
mula 

Pb 1 _ 2 .La z [(Mg^Nb { _ATi 1 _*] l-z/4^3+x/2-3xt/2-xz/8+3xtz/8 

wherein x is between about 0.6 and about 0. 95 , z is between 
about 0 and about 0.08, and t is between about 0.30 and about 
0.36. In especially preferred electro-optic ceramic materials 
of the present invention, z is between about 0.01 and about 
0.06. In preferred electro-optic materials of the present inven- 
tion, t may be between about 0.32 and about 0.34, and x may 
be between about 0.65 and about 0.90. 

According to the present invention, the preferred OptoCe- 
ramic® PMN-PT is lanthanum modified solid solution of 
lead manganese niobate and lead tantalate (La: PMN-PT, or 
PLMNT) with a nominal 3.5/75/25 La/PMN/PT composi- 
tion. For reference purposes nominal values for this PLMNT 
are n=2.48 and R=1 1 .5xl0 -16 m 2 /V 2 at 1 .06 jam and n=2.458 
and R=6.86xl0 -16 m 2 /V 2 at 1 .55 pm, both at 23° C. PLMNT 
has a polycrystalline structure with crystal sizes ranging from 
about 5 to 20 microns. In the form of hot-pressed ceramics, it 
is optically isotropic and exhibits little birefringence with 
zero applied voltage. 

An electrode geometry of phase retarder which takes 
advantage of these materials is illustrated by the transverse 
field configuration illustrated in FIG. 3 . A plate 30 , for use as 
the variable phase retarder 12 and 13 of FIG. 1 , is shown 
which has thin metallized electrodes 31 and 32 on the both 
sides, respectively, of a block section 33 of EO material. The 
electrodes 31 and 32 are shown connected to a driving voltage 
34 provide by controller 102 shown in FIG. 1 . The electric 
field (designated E) is 90° to the direction of light propaga- 
tion. Since the effect is based on the electric field in the 
material, it is desirable to arrange the electrodes as close 
together as possible to minimize the control voltage. The 
optical axis of the EO element is perpendicular to the electric 
field. 

As an example, an EO phase retarder is made from a 
PLMNT plate with metallized electrodes, which has dimen- 
sions where the width (w) between the electrodes is about 5 
mm, a thickness (t) of about 4 mm, and a length (1) of about 5 
mm. Such a dimensioned PLMNT plate has an insertion loss 
of about 0. 1 dB and half-wave voltage of 320V for a 700 nm 
wavelength laser at 23° C. 

It will be understood that the particular geometry described 
is an example and that other device geometries can be used. 
Different dimensions of w and t can be selected to provide the 
same jt phase shift with different control voltage. The other 
preferred PLMNT plate has dimensions of, but not limited to 
a width (w) of about 0.5-5 mm, a thickness (t) of about 1 .30-5 
mm, and a length (1) of about 1 .0-5.0 mm. 

Another preferred OptoCeramic® material is PLZT with a 
nominal 9.5/65/35 La/Pb/Zr composition. This composition 
is known to be transparent in a range from 450 nm to 7 pm; 
see, for example, Lionel M. Levinson, Electronic Ceramics, 
Chapter 7 (Marcel Dekker, New York, 1987). 

Yet another preferred OptoCeramic® material, PZN-PT, 
for use in the present invention may be described by the 
general formula 
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P b 1 - y - z Ba > ,La z [(Zn / Nb { _^) x Ti x _ x \ { _ z / 40 3+x/2 -3xt/2-xz/s + 

3 xtz/S 

wherein x is between about 0.5 and about 0.9, y is between 
about 0.05 and about 0.5, z is between about 0 and about 0.05, 
and t is between about 0.30 and about 0.36. In especially 
preferred electro -optic ceramic materials of the present 
invention, x is between about 0.65 and about 0.85, y is 
between about 0. 1 and about 0.2, z is between about 0.02 and 

0.04, and t is between about 0.32 and about 0.34. 10 

Referring to FIG. 4 is the measured variable phase retarder 
response. The optical response 41 to the applied control volt- 
age 42 is less than 1 microsecond. This means that the speed 
of setting and resetting of the VPR is as fast as 1 .0 microsec- 
ond. According to the present invention, the imaging appara- 15 
tus only needs 2 and 3 settings for the first and the second 
variable phase retarders 12 and 13 to complete a measurement 
cycle, respectively. The polarization sampling time could be 
less than 3 microseconds. This equals the system can generate 
a set of 4 Stokes component images at a speed of up to 300K 20 
frames per second (fps). With a fast enough imaging sensor, a 
real-time Stokes polarization imaging system can be realized. 

In the present invention, another preferred EO material for 
variable phase retarders 12 and 13 are made from liquid 
crystals for a polarization imaging system working between 25 
the optical wavelengths 350 nm to 2000 nm. Using the 
invented control algorithm illustrated in FIG. 2, and relation- 
ship between sequential intensity images and Stokes images 
listed in Table 1 , a fast polarization imaging system using LC 
variable retarders can also be configured. 30 

A polarization imaging apparatus shown in FIG. 1 is basi- 
cally a passive imaging system which using ambient light 
source to illuminate the sample. Referring to FIG. 5 is a 
schematic view of a preferred embodiment of an active polar- 
ization imaging apparatus in accordance with the present 35 
invention. It consists of an illumination unit 50 and an imag- 
ing sensing unit 500. The said illumination unit and imaging 
sensing unit are aligned in a reflection configuration. Both 
illumination unit 50 and imaging sensing unit 100 can be set 
to an angle related to sample 56 surface normal 57 . In addition 40 
to a light source 51, the illuminator may consist of parts such 
as a collimator 52, a polarizer or polarization modulator 53, 
and a fixed wavelength filter or tunable wavelength filter 54. 
Tunable optical band-pass filter was used to produce desired 
wavelength for the imaging, from visible to infrared (IR). The 45 
illuminator could produce a linear or arbitrary polarized light 
beam. 

A preferred light source 51 is a broadband light source, 
such as a tungsten-halogen lamp which provides a strong 5Q 
intensity over a broad spectrum from UV to NIR (300 nm to 
2 pm). Another preferred light source is a light emitting diode 
(LED) source or laser source. 

The scattering light 58 from sample 56 is detected by the 
imaging sensing unit 500. The optical lens set 11 in 100 55 
shown in FIG. 1 can be an optical objective lens 59 in FIG. 5 
with an infinite focal length, which makes the polarization 
imaging system a Stokes microscope. 

The first VPR 12 and second VPR 13 were controlled 
independently by a computer 102 through a controller unit 60 
101 which generates a sequential of voltages to control the 
phase retardations ofVPRs 12 and 13. A set of four intensity 
images, I 0 , I x , I 2 and I 3 of the sample were captured by imag- 
ing sensor 15 when the phase retardations ofVPRs 12 and 13 
were set at (0,0), (tc, 0), (jt,jt) and (jt/2,jt), respectively Then 65 
four Stokes components of a Stokes image, S 0 , S 1; S 2 and S 3 
can be calculated by computer 102 according to equation (1). 
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Referring to FIG. 6 is a schematic view of a polarization 
imaging apparatus shown in FIG. 5 configured in a transmis- 
sion configuration. 

It will be apparent to those skilled in the art that various 
modifications and variations can be made to the present 
invention without departing from the spirit and scope of the 
invention. Thus it is intended that the present invention cover 
the modifications and variations of this invention provided 
they come within the scope of the appended claims and their 
equivalents. 

What is claimed is: 

1. A polarization imaging apparatus for measuring Stokes 
images of a sample comprising: 

an optic lens set, 

a first variable phase retarder optically coupled to the optic 
lens set and oriented 22.5-degree angle with respect to 
the linear polarizer, the first variable phase retarder hav- 
ing a first controllable phase retardation; 

a second variable phase retarder optically coupled to the 
first variable phase retarder, the second variable phase 
retarder oriented at a 45 -degree angle with respect to the 
linear polarizer, the second variable phase retarder hav- 
ing a second controllable phase retardation; 

a linear polarizer optically coupled to the second variable 
phase retarder; 

an imaging sensor optically coupled to the polarizer for 
sensing the intensity images of the sample through the 
linear polarizer; 

an controller for varying the first controllable phase retar- 
dation of the first variable phase retarder and for varying 
the second controllable phase retardation of the second 
variable phase retarders; and 

a computer for controlling the controller, recording a 
sequential of four intensity images by varying the first 
and second variable phase retarders and converting the 
four sequential intensity images to four components of a 
Stokes image. 

2. The polarization imaging apparatus of claim 1 wherein 
the four components of a Stokes image are converted from 
four sequential intensity images when setting the first and 
second variable phase retarders at four different pair of phase 
retardations. 

3. The four components of a Stokes image in claim 2 are 
converted from four sequential intensity images when setting 
the first and second variable phase retarders having the phase 
retardations of (0,0), (jt,0), (jt,jt) and (jt/2,jt), using the for- 
mula: 


r So = 0.5 x (/ 0 + h) 

Si =0.5x(/ 0 -/i) 

‘ 5*2 = 0.5 x (/ 2 - So) 

As = 0.5 x (/ 3 - So) 

4. The polarization imaging apparatus of claim 1 wherein 
the four components of a Stokes image are converted from 
four sequential intensity images when setting the first and 
second solid-state ferroelectric variable phase retarders at 
four different pair of phase retardations. 

5. The first and second solid-state ferroelectric variable 
phase retarders of claim 4 each comprises one of the materials 
in the lanthanum modified PMN-PT (PMN-PT) system, or 
lanthanum modified PZN-PT (PZN-PT) system. 

6. The polarization imaging apparatus of claim 1 wherein 
the four components of a Stokes image are converted from 
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four sequential intensity images when setting the first and 
second liquid crystal variable phase retarders at four different 
pair of phase retardations. 

7. A method to measure Stokes images of a sample, the 
method comprising the steps of: 

orienting a first variable phase retarder at a 22. 5 -degree 
angle with respect to the linear polarizer; 
orienting a second variable phase retarder at a 45 -degree 
angle with respect to the linear polarizer; 
sensing a sequential of four intensity images with four sets 10 
of phase retardations while varying the first controllable 
phase retardation of the first variable phase retarder and 
varying the second controllable phase retardation of the 
second variable phase retarders; and 
converting a sequential of four intensity images to four 15 
components of a Stokes image. 

8. The method in claim 7, wherein the four sets of phase 
retardations of first and second variable phase retarders hav- 
ing the phase retardations of (0,0), (jt,0), and (jt/2,Jt). 
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9. The method in claim 7, wherein the four components of 
a Stokes image are converted from four sequential intensity 
images when setting the first and second variable phase 
retarders having the phase retardations of (0,0), (or,0), (jt,jt) 

5 and (jt/2,jt), using the formula: 

' So = 0.5 x (/ 0 + /i) 

Si =0.5x(/o-/i) 

' 5 2 =0.5x(/ 2 -/ 0 ) 

^3=0.5x(/ 3 -/o) 

10. The method in claim 7 further comprising the steps of: 
illuminating the sample with a light source; 
adjusting the light source polarization state; and 
changing the light source wavelength spectrum. 



